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ABSTRACT 

The aim of this work is the determination of the twist, writhe, and self magnetic hehcity of penumbral 
filaments located in an inner Sunspot penumbra. To this extent, we inverted data taken with the 
spectropolarimeter (SP) aboard Hinode with the SIR (Stokes Inversion based on Response function) 
code. For the construction of a 3D geometrical model we applied a genetic algorithm minimizing the 
divergence of B and the net magnetohydrodynamic force, consequently a force-free solution would be 
reached if possible. We estimated two proxies to the magnetic helicity frequently used in literature: 
the force- free parameter az and the current helicity term hc^. We show that both proxies are only 
qualitative indicators of the local twist as the magnetic field in the area under study significantly 
departures from a force-free configuration. The local twist shows significant values only at the borders 
of bright penumbral filaments with opposite signs on each side. These locations are precisely correlated 
to large electric currents. The average twist (and writhe) of penumbral structures is very small. The 
spines (dark filaments in the background) show a nearly zero writhe. The writhe per unit length of 
the intraspines diminishes with increasing length of the tube axes. Thus, the axes of tubes related to 
intraspines are less wrung when the tubes are more horizontal. As the writhe of the spines is very 
small, we can conclude that the writhe reaches only significant values when the tube includes the 
border of a intraspine. 

Subject headings: methods: observational - methods: numerical - Sun: magnetic topology - sunspots 
- techniques: polarimetric 
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1. INTRODUCTION 

Investigating the physical nature and the dynamics of 
penumbral filaments is essential in order to understand 
the structure and the evolution of sunspots and their sur- 
rounding moat regions. Many important observational 
aspects of penumbral filaments are well settled down, al- 
though their interpretation is still a source of debate, 
e.g., the brightness of penumbral filaments, the inward 
motion of bright penumbral grains, the Evershed flow, 
the Net Circular Polarization (NCP), as well as mov- 
ing magnetic features in the sunspot moat. During the 
last few years, new observational discoveries, e.g., dark 
cored penumbral filaments (Scharmer et al. 2002), strong 
downflow patches in the mid and outer penumbra (Ichi- 
moto et al. 2007a), penumbral micro-jets (Katsukawa et 
al. 2007; Jurcak & Katsukawa 2008), or twisting motions 
of penumbral filaments (Scharmer et al. 2002; Rimmele 
& Marino 2006; Ichimoto et al. 2007b; Ning et al. 2009) 
broadened the number of unknowns and gave new im- 
pulse to the investigation of sunspots. For recent reviews 
see Borrero (2011), Bellot Rubio (2010), Borrero (2009), 
Schlichenmaier (2009), or Tritschler (2009). An espe- 
cially controversial issue is the study of the twisting mo- 
tions of penumbral filaments. On the one hand, Ichimoto 
et al. (2007b) consider twisting motions as an apparent 
phenomenon, produced by lateral motions of intensity 
fluctuations associated with overturning convection. On 
the other hand, Ryutova et al. (2008) propose that the 
observed twist is an intrinsic property of penumbral fila- 
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ments and is produced as a consequence of rcconncction 
processes which take place in the penumbra. Su et al. 
(2008, 2010) conclude that the twist of penumbral fila- 
ments changes with time caused by an unwinding pro- 
cess. 

In the present paper we study the twist of filaments of 
the inner penumbra of a sunspot by means of the mag- 
netic helicity. We take advantage of the 3D geometrical 
model of a section of the inner penumbra of a sunspot de- 
scribed in Puschmann et al. (2010a) [hereafter. Paper I]. 
We use observations of the active region AR 10953 near 
solar disk center obtained on 1*** of May 2007 with the 
Hinode/SP. The inner, center side, penumbral area un- 
der study was located at an heliocentric angle 6' = 4.63°. 
To derive the physical parameters of the solar atmo- 
sphere as a function of continuum optical depth, the SIR 
(Stokes Inversion based on Response function) inversion 
code (Ruiz Cobo & del Toro Iniesta 1992) was applied to 
the data set. The 3D geometrical model was derived by 
means of a genetic algorithm that minimized the diver- 
gence of the magnetic field vector and the deviations from 
static equilibrium considering pressure gradients, gravity 
and the Lorentz force. We can not assess the unicity of 
the resulting model: the found solution just minimizes 
the divergency of the magnetic field and the modulus of 
the net force, neglecting the contribution of the aceler- 
ation terms. For a detailed description we refer to Pa- 
per I. In Puschmann et al. (2010b) [hereafter. Paper II], 

we calculated the electrical current density vector J in 
the above mentioned area and found the horizontal com- 
ponent of the electrical currents 4 times larger than 
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Fig. 1. — Left panel: Length i of the field lines of B integrated 
from the top (z = 200 km) to the bottom layer (z = km) or until 
X = 4.2 Mm. Thick isolines of equal I define 14 areas (flux tubes) 
related to spines (black) and intraspines (white), respectively. Thin 
isolines denote different thresholds for selecting the flux tube area. 
Right panel: Vertical component oi B at z = 200 km. Contour lines 
correspond to horizontal cuts through the 14 flux tubes with the 
largest area at 2 = 200,150,100,50,0 km (the thickness of each 
line diminishes with depth). Black and white contours distinguish 
between flux tubes related to spines and intraspines, respectively. 



the vertical component (thus confirming the results of 
Pevtsov & Peregud 1990; Georgouhs & LaBonte 2004). 
In addition, we concluded that the magnetic field at the 
borders of bright penumbral filaments departs from a 
force-free configuration (see also Zhang 2010). These re- 
sults are strongly significant considering that we have 
imposed that our solution minimizes the net force, in- 
cluding the Lorentz force, and consequently, a force free 
solution should be found if it were possible. 

We can evaluate the magnetic field lines by the inte- 
gration of B starting at each pixel of the top layer of 
our volume of the inner penumbra (of 4.2 Mm x 5.6 Mm 
X 0.2 Mm, see Paper II). In the left panel of Fig. 1 wc 
present the length £ of each field line. As our sunspot 
has negative polarity, B points downwards and wc thus 
integrate the field lines from the top layer to the bottom 
layer. The majority of field lines end up at the bottom 
layer, except the lines starting at larger X coordinates 
in our FOV. In the whole analized volume the magnetic 
field has the same polarity, and consequently, the field 
lines do not present maxima nor minima in our region, 
i.e., the field lines always travel downwards. This fact, as 
we will see later, simplifies the evaluation of the magnetic 
helicity. Areas with larger i correspond to intraspines, 
since B is more horizontal, the field lines thus traverse 
larger distances inside our volume. 

In the left panel of Fig. 1, we selected 14 areas, 7 cor- 
responding to intraspines (thick white contours) and 7 to 
spines (thick black contours), according to the length i 
of the magnetic field lines. For each of the selected zones 
we define a volume delimited by the field lines setting off 
from each pixel of the closed curve of the top layer. In 
the right panel of Fig. 1 we show the vertical component 



oi B a.t z = 200 km together with horizontal cuts through 
each of these volumes at z — 200, 150, 100, 50, & 0km. 
Since the umbra is placed at the right hand side in the 
FOV, the cuts at deeper layers are displaced to the right. 
Finally we checked that the magnetic flux traversing each 
of the cuts is approximately constant: the standard de- 
viation of the relative variation of the magnetic flux be- 
tween the top and the bottom layer is 4.5%. Thus, each 
volume can be approximately considered as a flux tube. 

However, the 14 areas were selected quit arbitrarily. 
In order to study the dependence of twist, writhe and 
magnetic helicity on the fiux tube area, 40 additionally 
smaller tubes have been defined inside the larger ones 
(thin isolines in the left panel of Fig. 1). Thus, for most 
of the 14 zones wc have several tubes of different size, 
many of them being the internal part of the larger one. 
Consequently, we have 54 tubes, 27 of them related to 
intraspines and 27 to spines. 

2. MAGNETIC HELICITY, TWIST, AND WRITHE 

The study of helicity of solar magnetic features has 
been a hot topic during at least the last 25 years. Mag- 
netic helicity has been investigated in solar structures at 
different spatial scales in the photosphere and chromo- 
sphere, as well as in the solar wind (see e.g. the reviews 
of Brown et al. 1999; Rust 2002; Pevtsov & Balasubra- 
maniam 2003; Dcmouhn 2007; Dcmoulin & Pariat 2009, 
and references therein). The helicity in penumbral fila- 
ments has been analyzed by means of some proxies by, 
e.g., Ryutova et al. (2008), Tiwari et al. (2009), Su ct al. 
(2010), and Zhang (2010). 

The magnetic helicity, H„i, quantifies how the mag- 
netic field is twisted, writhed, and linked. Hm plays a 
key role in magneto-hydrodynamics because it is almost 
conserved in a plasma having a high magnetic Reynolds 
number (see e.g., Berger 1984). The magnetic helicity 
of a vector field B, fully contained within a volume V 
and bounded by a surface S (i.e., the normal component 
Bn — B-n vanishes at any point of S), is (Elsasser 1956): 



A ■ B d-^x, 



(1) 



where the vector potential A satisfies B ~ V x A. Berger 
& Field (1984) showed that Eq. 1, is not gauge-invariant 
if the volume of interest is not bounded by a magnetic 
surface, i.e., if B crosses S (as in the case of the volume 
of the penumbra retrieved from our observations). In 
this case the relative magnetic helicity (Finn & Antonsen 
1985) should be used: 



H. 



rel 



(A + Ap) ■ (B - Bp) d^x, 



(2) 



where Bp is a potential field having the same normal 

component i?„ on 5, and Ap is its vector potential. 
The relative helicity reflects twist, writhe, and linkage 
with respect to a current-free (potential) fleld, i.e., its 
minimum-energy state for the given i3„-condition on 
S. The relative magnetic helicity so deflned is gauge- 
invariant and has the same conservation properties and 
amount of topological information as the magnetic helic- 
ity. Throughout this article, the term magnetic helicity 
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Fig. 2. — From left to right: local twist (T'""), az/in, current helicity density {he), and Ahc^ evaluated at z ■ 
the same as in the left panel of Fig. 1. 
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refers to the relative magnetic hehcity. For an isolated 
magnetic flux rope, H^f is proportional to the sum of 
its twist T and writhe W (Bcrgcr & Field 1984; Torok 
et al. 2010): 

H;:' = {T + (3) 

where $ is the magnetic flux of the rope. The twist 
is the turning angle of a bundle of magnetic field lines 
around its central axis, whereas the writhe quantifies the 
helical deformation of the axis itself. Following Berger & 
Prior (2006), the twist of an infinitesimal rope is given 
by T = J T'"^/, I being the arc length along the central 
field line of the rope and T^"" the local twist: 



dl inB 



II 



(4) 



being J| | and Bu the components of the current and mag- 
netic field parallel to the central field line of the rope. 
With this definition, T = 1 when the field lines twist 
around the axis by an angle of 2tt and T^°° is the local 
twist per unit length evaluated at a given geometrical 
height at each pixel. If the rope has a non-infinitesimal 
cross section E, the local twist of the rope Tl°° is given by 
the average of the infinitesimal local twist, T'"" over E. 
Berger & Prior (2006, sec also Torok et al. (2010)) give 
expressions for the writhe of specific geometrical con- 
figurations. Provided that the magnetic field lines in 
the inner penumbral region studied in this paper always 
travel downwards in z, i.e. without showing maxima nor 
minima between the zq =200 km and zi =Okm height 
layers, we can evaluate the writhe using a very simplified 
formula: 



1 r^i 
W= — 
27r 



1 



1 



-(f X r')^dz. 



(5) 



obtained as a particular case of the more general for- 
mula of Berger & Prior (2006). In Eq. 5, r stands for 



the tangent vector to the tube axis; for the vertical 
component of r ; and t' = ^ . 

Equations 3, 4, & 5 allow the evaluation of the twist, 
the writhe and the magnetic helicity of an isolated tube. 
What happens in the case of a non-isolated tube, as is 
clearly the case of the penumbral tubes? The writhe of 
a magnetic rope, isolated or not, is a measure of the he- 
lical deformation of the axis of the rope, while the twist 
quantifies the winding of the magnetic field lines of the 
rope around its axis. The magnetic helicity measures the 
linking number of the field lines, averaged over all pairs 
of lines, and weighted by the flux (Berger & Prior 2006; 
Moffatt 1969). We can simplify the case of a non-isolated 
tube to a scenario in which we have just two adjoining 
tubes. It is clear that we can define the writhe and twist 
for each individual tube, and consequently its magnetic 
helictiy, but the helicity of the whole configuration is not 
just the sum of both contributions. We would need to 
include an extra term taking into account the linking be- 
tween both tubes: the mutual helicity. Consequently the 
application of equations 3, 4, & 5 to our tubes retrieves 
only the contribution of the local values of the twist, 
writhe and the self helicity. The contribution of the sur- 
rounding tubes to the helicity of each flux tube is not 
considered. This contribution, the mutual helicity, could 
be larger than the self helicity (see e.g., Rcgnier & Priest 
2007). The mutual helicity can be calculated using the 
procedure described in Berger & Prior (2006). However, 
for the scope of this paper we limit the calculation to self 
helicity. 

3. PROXIES OF THE MAGNETIC HELICITY 
The current helicity density is defined (see e.g., See- 
hafer 1990) as he = B ■ V x B. If we take as mag- 
netic ropes the tubes defined by the field lines start- 
ing in each pixel, (being then = B), Eq. 4 becomes 
r'°° = hc/4:TTB^. The parameter a is usually defined in a 
force- free configuration, i.e., when B is parallel to its curl. 
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by V X _B = aB. The a parameter can be defined for non 

force-free fields: a = B-VxB/B^. Tliis is tlie definition^ 
we will use throughout the paper. In order to see how 
this re-defined a differs from the force-free definition, let 
us decompose V x B = (V x _B)|| + (V x B)^. From 

its definition, a becomes equal to a ~ ±|(V x B)^^\/B 
which is equal to the classical definition for a force-free 
field. The ± is needed to consider the case when ( V x B) 1 1 

and B point in opposite direction, i.e., when a is nega- 
tive. For a no force free-field, the parameter a is then the 
ratio between the parallel component of the curl of the 
magnetic field and its modulus. Evidently, in a general 
case, we will have a = hc/B^ = AnT^"". 
Given the difficulty of empirically obtaining and 

J, one finds many works where different proxies were 
used. Before evaluating the magnetic helicity we can 
calculate, from our data, some of the most usual proxies 
of the magnetic helicity. Among them the most common 
proxies, with several different but more or less equiva- 
lent definitions, are the az = {V x B)z/Bz parameter 
(see e.g., Su et al. 2009, 2010; Pevtsov et al. 2008, and 

references therein), and the parameter ~ i?z(V x B)^ 
(see e.g., Zhang 2010). It is evident that for a force-free 
field az = a and = aB^ = hcB^/B^. That means 
that the parameter could be meaningless for nearly 
horizontal magnetic fields, such as those found in sunspot 
penumbrae. 

Many authors suppose that the sign of the integral of 
Uz (or hc^ ) over the volume of a magnetic structure coin- 
cides with the sign of iJ^'Jf ' , although this fact has not yet 
been demonstrated (Dcmoulin 2007). On the other hand, 
Hagyard & Pevtsov (1999) point out that hc^ only con- 
siders the vertical component Jz of the electric current 
density vector, hc^ can strongly differ from he, provided 
that Jz is much smaller than the horizontal components, 
at least in the inner penumbra (see Paper II). Besides, 
Pariat et al. (2005) comment that, since the magnetic he- 
licity is a global quantity, it is not obvious that a helicity 
density has any physical meaning. 

We evaluated these proxies in the inner penumbral 
region under study. In the left panels of Fig. 2 we 
present T^"" (evaluated from Eq. 4) and az/47r evalu- 
ated at each pixel at z = 200 km. As in the inner 
penumbra the magnetic field is not force-free at the bor- 
ders of bright penumbral filaments (see Paper II), az/4:TT 
(2"'^ panel) is only qualitatively similar to T'°°. In the 

3'''^ and 4*^ panel we present he ^ B ■ V x B and he, 
evaluated at z = 200 km. he, is multiplied by a fac- 
tor 4 just to make easier its comparison with he- As 
we have seen before, T'°" = he/AnB^, and thus the 
general aspect of he is very similar to T"'"". However, 
he^ resembles he only marginally (the standard devia- 
tions are a{he) = 2.1 G^m^^ and <T{heJ = OAG^ in~^). 
The values of az/4:TT and he, at z = 200 km obtained 
here are very similar to the results found in the lit- 
erature: Tiwari et al. (2009) found that az/47r varies 
around ±0.15Mm~^ along azimutlial paths in the mid- 
dle penumbra; Su et al. (2010) found a fluctuation of 

^ Yeates et al. (2008) denominate current helicity to this gener- 
alized a parameter. 



Uz/'Itt larger than ±0.05Mto~^ over an inner penumbral 
region; Su et al. (2009) found that he^ fluctuates along an 
azimuthal path in the inner penumbra with an amplitude 
larger than 1 while Balthasar & Gomory (2008) 

found penumbral mean values of about 0.04G^to~^. In 
the four panels of Fig. 2, the outlined areas were selected 
by different thresholds of £, the length of the magnetic 
field lines between the layers z = 200 km and z = Okm 
(see also Section 1 and Fig. 1). The sign of the inte- 
gral of T'"" and az/4:T: over the above mentioned areas 
only coincides in 39% of the 54 tubes (if we consider only 
the areas related to the intraspines this value decreases 
to 15% of the 27 tubes). The same figures are obtained 
for the percentage of coincidence between the signs of the 
integrals of he and he^ , approximately. This weak coinci- 
dence demonstrates that, at least in the inner penumbra 
of a sunspot, az/47r and hc^ are not good estimates of 
T'"" and he, respectively. As already shown in Paper II, 
the magnetic field in the area under study significantly 
departures from a force-free configuration. 

Note that T'"" reaches significant values only at the 
borders of the intraspines: these are exactly the areas 
where the electric current density is large (see Fig. 1 
of Paper II). Furthermore, often T'°° changes its sign at 
both sides of bright filaments, i.e., for the majority of 
the intraspines, T'"" shows negative values at the upper 
(larger Y-coordinate) borders of the filaments and posi- 
tive values at the lower borders. The alternation of the 
sign in the twist is also observed (although less evident) 
in the maps of az and he^ and it is clearly visible in Ti- 
wari et al. (2009), Su et al (2009), Su et al. (2010), and 
Zhang (2010). 

This phenomenon can be explained if we consider that 
the field lines of the magnetic background component 
wrap around the intraspines (Borrero et al. 2008) and 
tend to meet above the intraspines, thus generating a 
curvature of different sign in the field lines at both sides 
of the intraspines. Thus T'"" could be measuring the 
twist of the background field wrapping around the in- 
traspines rather than the twist of the field lines of the 
intraspines themselves. However, the alternation of signs 
of the twist at both sides of a penumbral filament is com- 
patible with the magnetohydrostatic equilibrium model 
of a magnetic fiux tube built by Borrero (2007). This 
model includes a transverse component of B having op- 
posite twist at both sides of a plane longitudinally cutting 
the fiux tube. This model is able of explaining both the 
dark cored penumbral filaments and the net circular po- 
larization observed in penumbral filaments (Borrero et 
al. 2007). Magara (2010) suggests the existence of an 
intermediate region where the magnetic field has a tran- 
sitional configuration between a penumbral fiux tube and 
the background field: in such areas, coinciding with the 
largest electrical current density (see Paper II), penum- 
bral micro-jets are produced as observed by Katsukawa 
et al. (2007). 

4. NUMERICAL TEST 

To check its correctness, the procedure used for the 
evaluation of twist, writhe, and magnetic helicity was 
applied to two different analytical cases. In the first case 
we consider that the magnetic field lines follow a helix 
around a vertical straight line. The magnetic field vec- 
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tor is defined hy B = Bq z + Bir 9 with Bq and Bi 
being constant. B could easily be decomposed in a po- 
tential Bp = Bq z and a close (toroidal) field Be ^ Bir 9. 

Obviously, the potential component Bp fulfills the con- 
ditions required in Eq. 2 in the case of a cylindrical 
tube: B and Bp have the same normal component on 
the external surface of the tube. The determination 
of the vector potential, in this case, is straightforward: 
A = Bor/2 9 - Bir^/2z. The vector potential of the 

potential component will be Ap = Bq r/2 9. Using Eq. 2, 
the magnetic helicity of a cylindrical tube of height L and 
radius R becomes Hm = Bq Bi L. The magnetic 
flux of this tube is $ = ttBqR^. This tube has a zero 
writhe because its axis is a straight line. From Eq. 3, the 
twist follows as 

T=^. (6) 

This is obviously the expected result, provided that the 
pitch (of screw-step) of our helix is and the twist is 

a measure of the number of turns done by the magnetic 
field lines along a longitude L. 

In the first four rows of Table 1 we present the analyti- 
cal (i.e., using Eq. 6) and numerical results (using Eqs. 4 
and 5) for tubes with Bq = -0.2 T, Bi = 0.04TMm"\ 
L ~ 0.225 Mm and a radius R equal to 0.2 and 1 Mm. 
We used the same spatial grid as in the observational 
case. 

In order to check the accuracy of the determination 
of twist, writhe and magnetic helicity in a more gen- 
eral case, we carried out a second test, building a heli- 
cal magnetic field that turns around a helical axis. Let 
us suppose that the axis of the tube is a vertical he- 
lix of radius Rh turning an angle ^ through a length 
L = 0.22 Mm. Then, the magnetic field at the axis will 
he B = Bqz + BiRh9 with Bi = Bq -^/L. Following 
Berger & Prior (2006) , the writhe of a magnetic fiux tube 
whose axis is a helix can be easily evaluated in terms of 
its polar writhe (i.e., area/27r of the section of the unity 
sphere limited by the tantrix curve and the north pole. 
The tantrix curve is the path, the tip of the tangent vec- 
tor takes on the unit sphere). In our case, the writhe 
becomes: 

(7) 



2lT 



\B\ 



Once we have the axis, we can easily build a tube with 
a given twist around such an axis. We chose the radius 
of the tube as R — 0.2 Mm, and the radius of the helical 
axis as Rh = 0.4 Mm, and Bi = 0.5 TMm-^ The angle 
^ takes a value of 9.124 degrees in order to have an 
analytical writhe (i.e., using Eq. 7) of 0.001. The added 
twist takes the values 0.0, 0.001, 0.01, 0.1, and -0.001. 
Note that, as the magnetic helicity is proportional to the 
sum of twist and writhe, in the last case we will have 
a null magnetic helicity. The results of these tests are 
presented in the ten bottom rows of Table 1 . 

5. RESULTS 

In Table 2 we present the resulting values of the axis 
length, magnetic fiux, twist, writhe and self magnetic he- 
licity for the flux tubes of the 14 largest zones. As the 



TABLE 1 

Test results: Writhe (W), twist (T), and magnetic 
helicity {hm), for a helical magnetic field (first four 

rows) and a helical MAGNETIC FIELD WINDING AROUND A 
HELICAL AXIS. 



R [Mm] 
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Hm [Mx^] 


0.2 


analytical 


0.00 


-7.16e-3 


-4.52e-|-34 




numerical 


-9.e-28 


-7.16e-3 


-4.50e-|-34 
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analytical 


0.00 
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numerical 
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analytical 
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Fig. 3. — Panel (a): Length of the axis of each flux tube. The 
values corresponding to tubes of the same zone are connected by 
lines. Cross symbols correspond to intraspines and small circles to 
spines. Panels (b), (c), and (d): writhe, twist, sum of writhe and 
twist. 



magnetic helicity depends on the square of the magnetic 
flux, and our selected areas are very different in area, 
the resulting magnetic helicity varies over a wide range 
of several orders of magnitude. To study the dependence 
of the precedent quantities on the length of the respec- 
tive axis, in Fig. 3, we plot the length of the axis of each 
flux tube C, the writhe, the twist, and the sum of twist 
and writhe for the 54 selected tubes related to intraspines 
(index ranging form to 26) and spines (index ranging 
form 27 to 53). The values corresponding to tubes of 
the same zone (see left panel of Fig. 1) are connected by 
straight lines. As the magnetic fleld in the spines is more 
vertical than in the intraspines, the length of the axis of 
the tubes related to the spines is clearly shorter. For each 
intraspine/spine zone the length of the axis grows with 
the index because each of the related tubes was chosen 
in the interior of the preceding one. The writhe of the in- 
traspines does not follow a clear pattern, but most of the 
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Fig. 4. — Panels (a) and (b): normalized twist (T/C,) and writhe 
{W/Q versus the length of the axis of each fiux tube C^. The ab- 
solute value has been overplotted in panel (b) with triangle 
symbols. The dashed line is the linear fit of |Vy/CI- Panel (c): 
wavenumber. Panel (d): T/C, versus a/ATv averaged over the sec- 
tion of each flux tube at 2 = 200 km. Panels (e) and (f): T versus 
the average of a^/i-K and h^^ respectively. In panels (d), (e) and 
(f) the straight line with slope 1 has been overplotted. Cross sym- 
bols correspond to intraspines and small circles to spines. 



TABLE 2 

Axis length magnetic flux $, twist, writhe and 

MAGNETIC HELICITY FOR THE 14 LARGEST ZONES. ThE 7 FIRST 
(last) ROWS ARE RELATED TO INTRASPINES (SPINES). 



index 


C [Mm] 


4" [Mx] 


T 


W 


Hm [Mx2] 





0.79 


-3.05C-I-18 


-0.0138 


0.0258 


l.lle-l-35 


4 


0.80 


-1.56C-I-18 


-0.0272 


0.0182 


-2.19e-|-34 


9 


1.06 


-4.71e-|-18 


-0.0210 


0.0024 


-4.12e-|-35 


15 


0.92 


-2.55C-I-18 


-0.0224 


-0.0163 


-2.51e-|-35 


20 


0.80 


-5.690-1-18 


0.0111 


0.0041 


4.94e-|-35 


23 


1.02 


-2.640-1-17 


-0.0054 


0.0166 


7.80e-|-32 


25 


1.00 


-5.13C-I-17 


0.0096 


0.0079 


4.59e-|-33 


28 


0.32 


-4.43c-(-18 


-0.0023 


0.0000 


-4.42e-^34 


29 


0.33 


-3.410-1-18 


-0.0028 


0.0021 


-8.34e+33 


30 


0.34 


-6.290-1-18 


-0.0041 


0.0011 


-1.20e+35 


36 


0.34 


-7.730-1-18 


-0.0065 


0.0010 


-3.30e+35 


40 


0.37 


-6.30C-I-18 


-0.0041 


0.0015 


-l.Ole+35 


41 


0.36 


-1.13C-I-19 


-0.0001 


0.0007 


7.56e-|-34 


46 


0.35 


-2.04C-I-19 


-0.0010 


-0.0017 


-1.15e+36 



tubes have a positive writhe. The spines show a nearly 
zero writhe. The twist of the intraspines, however, takes 
nearly always negative values. Consequently, the twist is 
partially canceled by the writhe, thus the absolute value 
of the self magnetic helicity is, nearly always, lower than 
the absolute value of the twist. 
In panels (a), and (b) of Fig. 4 we plot the twist, and 



writhe per unit length versus the length of the axis of 
each flux tube C,. The normalized twist does not show a 
clear dependence with the length of the axis, but the ab- 
solute value of the normalized writhe in the intraspines 
clearly decreases with increasing The axes of tubes 
related to intraspines are less wrung when the tubes 
are more horizontal (i.e., in the central part of the in- 
traspines). As the writhe of the spines is very small, we 
can conclude that the writhe reaches only significant val- 
ues when the tube includes the border of a intraspine. 
In panel (c) we plot the wavenumber [1/ pitch), i.e., the 
number of turns done by the magnetic field per length 
unit as a function of Q. Flux tubes related to spines show 
slightly smaller values but any clear dependence is not 
observed. 

In panel (d) we plot the normalized twist versus a/47r 
(i.e., the local twist T^"" at each pixel) evaluated at 
z ~ 200 km and averaged over each structure. Given 
the good correlation between both magnitudes we can 
use the average of the local twist as a good proxy for 
the normalized twist. Consequently, we can explain the 
small obtained twist values in terms of the local twist: as 
we have seen in Fig. 2, T'"" reveals significant values with 
opposite sign at the borders of the pcnumbral filaments, 
leading to a cancellation of the twist when integrating 
over the filament. 

In order to assess the reliability of the most used prox- 
ies we plot, in panels (e) and (f), the average twist of each 
structure versus the average value of az/4:T: and hc^, re- 
spectively. Both proxies are very bad indicators of the 
average twist of intraspines but they give a qualitatively 
good agreement (better in the case of he J for the spines. 
This asymmetry could be explained by the fact that both 
az and hc^ are only related with the vertical component 
Jz of the electric current density vector, which, as we 
show in Paper II, is much smaller than the horizontal 
components, mainly in and around the intraspines. On 
the other hand, in a force-free configuration = a and 
he, = he, and then, the discrepancy between these pa- 
rameters is a clear result of the non- validity of the force- 
free approximation in the inner penumbra: In Paper II 
we have already shown that, at the borders of bright 
penumbral filaments, the magnetic field strongly depar- 
tures from a force-free configuration. 

6. CONCLUSIONS 

In the present work we calculated the parameter a and 
its proxy az , the current helicity density he and its proxy 
he, the twist, the writhe, and the magnetic helicity of 
different structures of the inner penumbra of a sunspot. 
The parameters are evaluated from a three-dimensional 
geometrical model obtained after the application of a ge- 
netic algorithm on inversions of spectropolarimctric data 
observed with Hinode (sec Puschmann ct al. 2010a, Pa- 
per I). We demonstrate, that in the inner penumbra the 
frequently used proxies az and hc^ are only qualitative 
indicators of the local twist (twist per unit length, evalu- 
ated under the assumption that the axis of a fiux tube is 
parallel to the magnetic field) of penumbral structures. 
As shown in (Puschmann et al. 2010b, Paper II), the 
magnetic field in the area under study many times de- 
parts significantly from a force- free configuration and the 
horizontal component of the electrical current density is 
significantly larger than the vertical one. 
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The local twist shows only significant values at the bor- 
ders of bright penumbral filaments and reveals opposite 
sign at each side of the bright filaments. The opposite 
sign might be the reason for a cancellation of the twist 
when integrating over the filament, thus the twist of the 
penumbral structures is very small. Significant values 
of the local twist are exactly related to areas where the 
electric current density is large. The local twist could 
be measuring the twist of the background field wrapping 
around the intraspines and/or the twist of the field lines 
of the intraspines themselves; in the latter case the in- 
ternal structure of the tube would consist in two " cotyle- 
dons" (at both sides of a vertical plane longitudinally 
cutting the tube), harboring each one a magnetic field 
of opposite twist, compatible with the MHS model of 
Borrero (2007). 

The writhe per unit length diminishes with increasing 
length (decreasing inclination) of the axis of flux tubes 
related the intraspines. The small amount of twist and 
writhe shown by the spines indicates that the background 
field lines, in these zones, are nearly straight. 

A future study should clarify if the helicity apparent 
in the intensity maps of penumbral filaments in the mid 
and outer penumbra of sunspots is produced by helical 
flux tubes with a strong writhe or just by spurious effects 
produced by lateral intensity fluctuations. In any case, 
it is clear that a twisted tube does not per se generate 
any intensity fluctuation similar to the observations by 



Ryutova et al. (2008). Rather, there is the necessity of 
a writhe of the tube, in such a way that different lon- 
gitudinal portions of the tube were at different optical 
depths producing changes in the observed intensity. We 
will extend the present work (and necessarily the work 
presented in Paper I and II) on the entire sunspot. 
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